The quenching of compressible edge states around antidots 
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We provide a systematic quantitative description of the edge state structure around a quantum an- 
tidot in the integer quantum Hall regime. The calculations for spinless electrons within the Hartree 
approximation reveal that the widely used Chklovskii et al. electrostatic description greatly over- 
estimates the widths of the compressible strips; the difference between these approaches diminishes 
as the size of the antidot increases. By including spin effects within density functional theory in the 
local spin-density approximation, we demonstrate that the exchange interaction can suppress the 
formation of compressible strips and lead to a spatial separation between the spin-up and spin-down 
states. As the magnetic field increases, the outermost compressible strip, related to spin-down states 
starts to form. However, in striking contrast to quantum wires, the innermost compressible strip 
(due to spin-up states) never develops for antidots. 

PACS numbers: 73.21.Hb, 73.43.-f, 73.23. Ad 



A quantum antidot is a potential hill in a two- 
dimensional electron gas (2DEG) usually defined by 
means of an electrostatic split gate, see Fig. In a per- 
pendicular magnetic field, electrons are trapped around 
the antidot in bound states formed by magnetic con- 
finement. Experimental studies of magnetotransport in 
quantum antidots reported over the last decade reveal a 
rich magneto-conductance structure in the quantum Hall 
regime Bllllllli,|l[li,|i[lS[I3 Some of the ob- 
served magneto-conductance features can be understood 
within a one-electron picture in terms of semi-classical 
and quantum electron dynamics [J. However, a majority 
of experiments confirm a central role played by electron 
interactions and spin effects in antidot measurements. 
This includes, for example, the first direct observation of 
the fractionally quantized electron charge fractional 
statistics 0| , the striking effect of the frequency doubling 
of the Aharonov-Bohm (AB) oscillations iSjlSl], the detec- 
tion of the Coulomb charging the observation of the 
Kondo effect ^3 f^nd selective spin-injection ^l|- Inter- 
est in antidot structures is also motivated by their poten- 
tial for spintronic applications, where they can be used 
to inject or detect spin-polarized currents |l2l | or even as 
quantum gates [l^ . The antidots also provide a system 
for investigating edge states in general, because of the de- 
tailed information one can obtain from the dependence 
of AB peak positions on field and gate voltage. 

A detailed microscopic understanding of antidot sys- 
tems therefore requires a rigorous theory accounting for 
both interaction and spin effects. In contrast to quantum 
dots and wires which have been the subject of intense the- 
oretical study (see e.g. ^3,^3 ^^"^ references therein), 
the energetics of quantum antidots has received, with 
only a few exceptions |l6l |. practically no attention. In 
particular, the central issue concerning the structure of 
edge states and the formation of the compressible strips 
around antidots still remains an open question. The 




FIG. 1: Calculated electrostatic confinement potentials 
Vconf{r) for a quantum antidot with different gate radii R. 
The inset shows a schematical layout of the device defined in 
GaAs heterostructure with a cap layer of the width c = 14 
nm, a donor layer of the width d = 36 nm, and a spacer layer 
of the width s = 10 nm. The donor concentration Ud = 6-10^'^ 
m~^. The above parameters (corresponding to the bulk elec- 
tron density ni,„ife ~ 2 • lO^^m"^) are used for all antidot 
structures discussed in the paper. 

structure of the antidot edge states represents an im- 
portant key to an understanding of various effects such 
as AB oscillations 0, , Coulomb charging Q and spin 
selectivity |l l| , and it has been a subject of recent lively 
discussions |j, . The main goal of the present paper 
is to provide a rigorous theoretical description for the 
spin-resolved structure of edge states around a quantum 
antidot. 

We consider an antidot defined within a GaAs het- 
erostructure by a circular gate of radius R as illustrated 
in the inset to Fig. ^ We assume that electron motion is 
confined to the plane parallel to the heterostructure inter- 
face, r = (x,y). The external electrostatic confinement, 
Vconfir) = Vschottky + Vd + Vg{r), includes the Schottky 
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barrier Vschottky = 0.8 eV and the potential due to a 
layer of donors of width d situated at a distance c from 
the surface 0|, Vd = — Ji^d (c + c?/2) e^/£o£r, with rid 
being the donor concentration and Sr the GaAs dielec- 
tric constant. The electrostatic potential Vg{r) due to a 
circular gate is given by an analytical expression provided 
by Davies (Eq. (3.17) in Ref. Il9|) . The electrostatic con- 
finement Vconf{T) for different gate radii R is shown in 
Fig. CI. 

Utilizing the circular symmetry of the structure we in- 
troduce cylindrical coordinates and write down the wave 
function in the form -i/'(r) = e*''^0(r), where I — 0, ±1, . . . 
is the orbital quantum number. The Schrodinger equa- 
tion in a perpendicular magnetic field B reads 
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where the lengths are measured in units of the lattice 
constant a (used for numerical discretization), energies 
in units of I2m* a}\ q = eBa^/h, and a ~ ±i de- 
scribes spin- up and spin-down states, We include 
electron interactions and spin effects within the frame- 
work of density functional theory (DFT) in the local 
spin density approximation (LSDA) The choice of 
DFT-I-LSDA for the description of many-electron effects 
is motivated, on one hand, by its practical implementa- 
tion efficiency within a standard Kohn-Sham formalism 
[2l| . and on the other hand, by the excellent agreement 
between the DFT+LSDA and exact diagonalization 
and variational Monte-Carlo calculations |Q, |2^ per- 
formed for few-electron systems. Within the framework 
of the DFT+LSDA, the total confinement potential can 
be written in the form 

Vir) = V^onfir) + Vnir) + V^,{r) + gfisBa, (2) 

where 
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is the Hartree potential including the contribution 
from mirror charges (6 is the distance from the 
2DEG to the surface), n{r) = "-"(''); ""l^) = 
J2i fFDiE—Ep) is the electron density, and/_FD 

is the Fermi-Dirac distribution. For the exchange and 
correlation potential V^^{r) we utilize a widely used pa- 
rameterization from Tanatar and Cerperly [2^ (see Ref. 
for the expHcit expressions for Vxcir))- This param- 
eterization is valid for magnetic fields corresponding to 
filling factors v > \, which sets the limit for the applica- 
bility of our results. The last term in Eq. Q accounts 
for the Zeeman energy where /ib = is the Bohr mag- 
netron, and the bulk g factor of GaAs is g = —0.44. We 



solve Eq. 1^ sclf-consistently expanding the wavefunc- 
tions into sin-basis. Because the antidot represents an 
open system, we choose the computational domain suf- 
ficiently large to ensure that the electron density in the 
bulk (i.e. far away from the antidot) is constant and does 
not change when we increase the domain size. 

Following our previous analysis of edge state structure 
in quantum wires liil , we start with the Hartree 
approximation [disregarding exchange and correlation in- 
teractions by setting V^^{r) = in Eq. lO ]. Figure|21(a) 
shows the electron density profiles (the local filling fac- 
tors) u{r) = n(r)/nB {ub = eB/h) around antidots with 
different radii for a representative value of magnetic field 
B = 4.2 T. The corresponding average wave function po- 
sitions ipi for different eigenenergies Ei (i.e. the magne- 
tosubbands) are shown in Figs. |21(b)-(d) illustrating the 
formation of the compressible strips around the antidots. 
[Note that in Fig. |2]each eigenstate V"; is represented by a 
dot. Because of a large number of the eigenstates, ^ 10^, 
the dots are merged into solid lines]. The compressible 
strips are composed of partially filled electron states that 
screen the external potential and lead to a flattening of 
the subbands in the compressible regions l^^. [Follow- 
ing Refs. 0, [2^ [2^ [23 We define the compressible 
strips within the window \E — Ep < 27rfcT|)]. Figure |21 
also shows the total confining potential V{r), Eq. ||2Jl. 
Note that the calculated spin-up and spin-down densi- 
ties and potentials are virtually indistinguishable on the 
scale of the figure. Hence, in this magnetic field interval 
the effect of the Zeeman term on the subband structure 
is negligible, so that we may refer to the Hartree results 
as being the case of spinless electrons. 

Figure 13 shows the width of the compressible strips 
for spinless electrons, it;^„p, around antidots of differ- 
ent radii. For a comparison, we also plot Chklovskii 



et al. 



j27| analytical expressions, 
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giving the 



length of the compressible strips at the edge of a semi- 
infinite 2DEG. [w^^p depends on two parameters, the 
depletion length Z, and the filling factor in the bulk, 
Vbuik = nbuik/n-B- We extract I from the calculated self- 
consistent density distribution by fitting to the depen- 

dence n(r) = ntuik { 7^ ) where Ubuik is the electron 
density far away from the antidot [2^ ] . It has 

been demonstrated that for the case of quantum wires the 
width of the compressible strips for spinless electrons cal- 
culated in the Hartree approach, is in very good 
agreement with the Chklovskii et al. predictions |27l | for 
Wcomp I2M Ell • This is obviously not the case for the an- 
tidot structures where the Chklovskii et al. predictions 
[irj greatly overestimates the compressible strip width. 
For example, for an antidot with the radius R = 200 nm, 
the innermost compressible strip (i.e. corresponding to 
the edge state closest to the antidot) starts to form at 
B K, 3.9 T, whereas for the 2DEG edge this strip starts 
to already form at B « 2.6 T, see Fig. O The difference 
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FIG. 2; (a) The electron density profile (local filling factor), 
v{r) — n{r)/nB for antidots of diflferent radii R calculated 
in the Hartree approximation, (b) The corresponding mag- 
netosubband structure (i.e. the overage position of the wave 
functions ipi). Fat solid lines indicate the total confining po- 
tential V{r), Eq. (P. Magnetic field B = 4.2 T. Temperature 

r = 1 K. 
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between w^^^p and w^^^p 



is most pronounced for small 
antidot radii R and decreases as R increases (note that 
the limit R ^ oo effectively corresponds to the case of 
a straight boundary, i.e. the semi- infinite 2DEG). The 
difference between w^^p and w'^^^ can be understood 
as follows. For the case of a semi-infinite gate an electron 
in the vicinity of the edge of the 2DEG experiences the 
Hartree potential originated from electrons in the semi- 
infinite region not covered by the gate. However, for the 
case of the antidot the Hartree potential is stronger as 
it includes an additional contribution from the electrons 
surrounding the antidot (that are otherwise depleted for 
the case of the semi-infinite 2DEG). This additional con- 
tribution effectively repels the electrons from the bound- 
ary towards the bulk of the 2DEG. This leads to less 
effective screening and thus to steeper potential prevent- 
ing the formation of compressible strips. 

Let us now analyze the spin-resolved edge state struc- 
ture within the DFT approximation. Figures 0] (c)-(f) 
show the electron density profiles (the local filling fac- 
tors) v'^ [r) = n'^{r)/nB and the average wave function 
position ■(/'f for different energies (the magnetosub- 
band structure) around the antidots with radius R — 200 
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FIG. 3: Width of the compressible strips around quantum an- 
tidot with different radii i? as a function of magnetic field cal- 
culated within the Hartree approximation (i.e. for the spinless 
electrons) a.nd its comparison to the Chklovskii el al. pre- 
dictions [23. A'^ = 1,2 and 3,4 refer to the subband number 
corresponding to the two lowest (spin degenerate) edge states. 
Temperature T = 1 K. 



nm for different representative magnetic fields B = 4.3 
T and B = 4.4 T. Within the Hartree approximation 
the subbands are virtually degenerate since the Zeeman 
splitting is very small in the magnetic field interval under 
investigation (see Fig. I^J. In contrast, the exchange in- 
teraction included within the DFT approximation causes 
the separation of the subbands for spin-up and spin-down 
electrons. Indeed, the exchange potential for spin-up 
electrons depends on the density of spin-down electrons 
and vice versa 0, 0| . In compressible regions the 
subbands are only partially filled (because fpD < ^ hi 
the the window \E — Ef\ < 2TTkT), and, therefore, the 
population of spin-up and spin-down subbands may be 
different. In the DFT calculation, this population differ- 
ence (triggered by Zeeman splitting) is strongly enhanced 
by the exchange interaction leading to different effective 
potentials for spin-up and spin-down electrons and even- 
tually to a subband spin splitting. As a result, the com- 
pressible region present in the Hartree approximation is 
suppressed and the spin-up and spin-down states become 
spatially separated by the distance dsep ~ 'W'rampi 
Figs. E](a)-(d). On further increasing the magnetic field 
the compressible strip starts to form for the outer (spin 



down) state such that d. 
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, see Figs. ^ 

(a)-(b),(e)-(f) {w'^^^p^ is the width of the compressible 
strip for the spin state a calculated in the DFT approx- 
imation). 

Far away from the antidot the subbands remain de- 
generate since they are situated below the Fermi energy 
E < Ep — 2TTkT and are thus fully occupied {fpD = !)■ 
As a result, the corresponding spin-up and spin-down 
densities are the same, hence the exchange and correla- 
tion potentials for the spin-up and spin-down electrons 
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FIG. 4: (Color online) (a) Width of the compressible strips for spinless electrons in the Hartree approximation, w^omp, as 
compared to dsep + Wc^mp' calculated within DFT. (b) With of the compressible strips for spin-up and spin-down states 
''^Imnjf I'lJ^cornif Calculated within DFT. (c), (e) The electron density profiles (local filling factor), v" {r) — n'^{y)/nB for different 
magnetic fields (indicated by arrows in (a)) B = 4.3 T and 4.4 T, and, (d), (f) the corresponding magnetosubband structure 
(i.e. the overage position of the wave functions ipf). Fat solid lines indicate the total confining potential V'^{r) , Eq. T = 1 
K. 



are equal, yj^(r) = Vj-^{r). 

A similar scenario for subband spin splitting also holds 
for quantum wires |2!tI |. However, an important and in- 
teresting distinction is that for quantum wires, as the 
magnetic field is increased, compressible strips form first 
for spin-up and then for spin-down states. In contrast, 
for the case of the antidot only the compressible strip 
for the spin-down state forms, whereas the compressible 
strip for the spin-up states (situated close to the antidot) 
never develops, see Fig. EJb) [This conclusion holds for 
all antidots sizes studied in this paper, sec Fig. Just 
as for the case of spinless electrons discussed above, we 
attribute this difference to less effective screening for the 
antidot structure leading to a rather steep potential near 
the antidot boundary that prevents formation of the com- 
pressible strip for the innermost (spin-up) state. Note 
that the absence of the compressible strip for the inner 
(spin-up) state is consistent with the interpretation of the 
unexpected doubling of the frequency of the Aharonov- 
Bohm oscillations explained in Ref. in terms of a 
charging of the outermost compressible regions around 
the antidots by electrons of the same spin. Our calcula- 
tions indicate that the charging of the innermost (spin-up 
states) is unlikely since the spin-up states do not form the 
compressible regions. [We note however, that wc arc not 
yet in a position to comment on whether charging of the 
compressible strips proposed in Ref. really does takes 
place. The answer to this question can be obtained from 
self-consistent trans por t calculations similar to those re- 
ported in e.g. Ref. |22|. Such calculations are currently 
in progress.] 

Finally we stress that all the results and conclusions 



presented in this paper for a temperature T = 1 K re- 
main valid for lower temperatures, since calculations per- 
formed for T = 0.2 K reveal that the width of the com- 
pressible strips remain practically unchanged. 

It is important to note that the effect of reduced screen- 
ing that strongly affects the antidot edge structure as 
discussed above might also be imperative for the case of 
quantum dot, where the spin selectivity in the edge state 
regime might strongly depend on the gate layout |3nl |. 
The implications of this effect for quantum dot geome- 
tries remains to be identified. 

To conclude, we find that for spinless electrons edge 
states around an antidot start to form for magnetic fields 
significantly higher than those predicted by Chklovskii et 
al. electrostatic description |27l |. By including spin ef- 
fects within spin density functional theory we show that 
the exchange interaction leads to qualitatively novel fea- 
tures in antidot edge state structure, such as the suppres- 
sion of compressible strips for lower fields, spatial sepa- 
ration between spin-up and spin down states, and to the 
total absence of the innermost compressible strip due to 
spin-up states (as opposed to the outermost compressible 
strip due to spin-down states that forms at higher fields). 
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